Dynamics of Selfavoiding Tethered Membranes II 
Inclusion of Hydrodynamic Interaction (Zimm Model) 
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The dynamical scaling properties of selfavoiding polymerized membranes with internal dimension 
D embedded into d dimensions are studied including hydrodynamical interactions. It is shown that 
the theory is renormalizable to all orders in perturbation theory and that the dynamical scaling 
exponent z is given by 2 = d. The crossover to the region, where the membrane is crumpled swollen 
but the hydrodynamic interaction irrelevant is discussed. The results apply as well to polymers 
{D = 1) as to membranes (D = 2). 



The statistical properties of polymerized flexible mem- 
branes, generalizing polymers, have found large interest 
during the last years. Due to the self avoidance, they are 
either found in a flat or crumpled swollen phase ^W ■ An 
analytical approach was initiated in pJs] , where calcula- 
tions of the static scaling exponent v were performed at 
1-loop order. Its consistency to all orders in perturbation 
theory has been established in jl^jl^. Recently, 2-loop 
calculations have been performed, which give reliable re- 
sults for all imbedding dimensions ||l^ , [l3| . 

In this letter we want to adress the question of the dy- 
namics of such membranes in the crumpled swollen phase 
embedded in some (viscous) solvent. (For a discussion of 
the flat phase, see 0.) 

We first summarize the main results before discussing 
the technical procedure to derive them. 

First of all, the Brownian motion of the particles, both 
of the solvent and of the fluid, have to be taken into 
account. It is responsible for the relaxation of the mem- 
brane. This can be studied via the auto-correlation func- 
tion which has for large membrane size and large time 
the scaling form 



D 1 




{{r{x,t)-r{x,0)f)^t 



2/z 



(1) 



Our goal is to determine z. 

In the physical system, hydrodynamics may be impor- 
tant. Two cases can be distinguished: In the first case, 
the hydrodynamic is irrelvant and the exponent z is given 
by 



+ df 



(2) 



where df is the fractal dimension of the membrane. This 
result has been established to all orders in perturbation 
theory in [[l5| . In the second case, hydrodynamic inter- 
actions are relevant and we will show below that this 
modifies the exponent z to 



= d 



(3) 



where d is the dimension of the embedding space, see be- 
low. This situation is plotted in figure 1, where also the 
phase-separation line is given. 



FIG. 1. The different phase-regions for a D-dimensional 
membrane embedded into d dimensions including hydrody- 
namique interactions. The region with 5 < and e < is the 
Gaussian phase. Selfavoidance and hydrodynamic interaction 
are irrelevant, i.e. v* — ^^J^ and z — j^- Hydrodynamics 
is relevant for small d and becomes irrelevant ii d > df + 2, 
where df = D/v* is the fractal dimension of the membrane. 
For d> df + 2 and e > 0, z = 2 + df. 

The exponent z should be observable via dynamic 
light- or neutron-scattering methods. To our knowledge, 
no such experiment has been performed. 

Theoretically, the dynamics for polymers has first been 
regarded in |lq,0] using scaling arguments. For mem- 
branes, a similar analysis has been performed in |Tq]. 

Scaling can best be studied using renormalization 
group methods. For polymers, such a treatment has been 
carried out at 1-loop order in |jlfl,0,n9^p4| . A proof of 
the renormalizability which ensures the correctnes of the 
method has been given in pq | for the case of purley dis- 
sipative motion (Rouse model). 

Let us now introduce the model with hydrodynamic 
interactions, show that it is renormalizable and calculate 
the scaling exponent z. 



The static behavior of the membrane is given by the 
Edwards Hamihonian 

nr] = \jj^^r{x))\hjj5'{r{x)~r{y)) . (4) 

The embedding of the Z)-diniensional membrane in d- 
dimensional bulk space is described by the mapping 
X g H ^- r(a;) G IR . & is the couphng constant, asso- 
ciated to the self-avoidance. 

Hydrodynamic interactions for polymers were first in- 
troduced by Zimm |25| . He wrote down the following 
Langevin equation, which we will also use for membranes: 



SiTT 

f{x,t)=V- (-_+C)(a:,i) 



(5) 



Here, • denotes the scalar product of the matrix operator 
V and the vector 6H/Sr, which is defined by 



(6) 



f -9-= / faix)ga{x) . 

The hydrodynamic interaction is 



'Dai3{x,y,r{x,t),r{y,t)) = 






(7) 

Jk{r{x,t)~r{y,t)) 



We will not repeat the derivation |£^ of equation (J^) 
here. Let us however note that one supposes that the hy- 
drodynamic degrees of freedom are fast enough, so that 
their dynmaics can be neglected and that screening ef- 
fects are irrelvant. This might be wrong for membranes 
and in this case our results would only apply to mem- 
branes with large holes. For 77 = 0, (H) reduces to purely 
diffusive motion (Rouse model). 
The noise correlation is 



Ca{x,t){V-C)0{y,t') = 25{t-t')S''{x~y) 



(8) 



This ensures that the static behavior is correctly repro- 
duced. Following Janssen |g6|, the corresponding field- 
theory is obtained by imposing the Langevin-equation 
through an auxiliary field f. Integration over the noise 
then yields the dynamic functional in Ito-discretization 



J = 



r ■ r + r ■ T> ■ — r ■ T> ■ r 

or 



This model has to be renormalized. Analogously to ITI] 
divergences only occur at small distances. They can 
be analyzed via a multilocal operator product expan- 
sion (MOPE). Renormahzability is ensured ||l|,|l^,|ll| if 
counter-terms for all possible marginal and relevant oper- 
ators are included into the action. Due to causality, only 
operators with at least one response field f are needed. 
These are the local operators ff^ f{—A)r and f^. Other 



local operators like rr" are forbidden by translation- 
invariance in r-space. 

By the same arguments one finds that there are no 
new marginal or relevant counter-terms proportional to 
2- or 3- body interactions and 4-body-interactions are ir- 
relvant. (Of course, long-range interactions are relevant, 
but they are not generated in perturbation theory.) We 
now want to show that the structure of the model is pre- 
served, i.e. that it can be renormalized if we introduce 
renormalized quantities according to 



A — ZxXr\I — 

r — y Zrn 
f = VZfji 



(10) 



,z„zz'z^/^-'^l' 



The two regularization parameters e and 6 are given by 



e^2D-vd 
5 = 2-vd 
2-D 



(11) 



V 



Perturbation theory is performed about the point ((5 = 0, 
e = 0), i.e. (Z) = 1, d = 4). As the model is constructed 
such that the static limit is correctly reproduced, Z and 
Zh are the renormalization-factors of the static theory 

Then, Z is determined in order to render r finite. 

The composition of f with any other operator with the 
same time argument is always free of divergences, as the 
contraction of f with any functional of r and f vanishes. 
We therefore conclude that the first term in the action, 
ff is correctly renormalized. 

Then, Z\ and Z^ are chosen to render the operator D 
finite. By the same arguments as above, the last term in 
the action, f ■ Vf is renormalized. 

We still have to show that also the composite operator 
I? • ^ is finite. This is a consequence of the equation of 
motion obtained through variation of J by f{x,t): 



f{x,t)+[V-^]{x,t) 



2{f-V){x,t)=0. (12) 



(9) This relation is valid as an operator-identity. We have al- 



ready renormalized r, D, -r— and f. Equation ( |12| ) thus 
states that also the composite operator I? • ^ is finite. 
We can now conclude that all the terms in the action are 
finite. 

There are three nontrivial relations which consider- 
ably simplify renormalization and which we are going 
to study now. First of all, due to the fact that in any 
interaction vertex the field r appears either as difference 
{r{x,t) — r{y,t)) or as spatial derivative, no divergence 



proportional to fr appears. (For a detailed discussion see 
flia.) This means that 



ZZ = 1 



(13) 



In addition, there is no proper renormalization of the hy- 
drodynaniic interaction. Let us explain this point. De- 
note by 



5al3 






^ik(r(x,t)~r(y,t)) ^ 



xfaix,t)gp{y,t) 



(14) 



any hydrodynamic interaction vertex. (The dotted line 
represents any polynomial in r and f or their deriva- 
tives.) Then singular configurations which give rise to 
a renormalization of the hydrodynamic interaction are 
those for which two interaction vertices are contracted to 
one single vertex. We claim that their multilocal opera- 
tor expansion (MOPE) does not contain a contribution 
proportional to the hydrodynamic interaction vertex: 



= 



(15) 



(The round dotted lines indicate points which are con- 
tracted.) This property is due to the analytic behavior 
of the long range (hydrodynamic) interaction. Dropping 
indices, the structure of such a contraction is 



We are interested in the IR-behavior. Suppose that the 
coupling related to the selfavoidance has flown to its IR- 
fixed point bji = 6*, what implies that also the scaling 
exponent v{bji), defined by 



vibn) 



2-D 1 d 
~^ 2'^dil 



InZ 



(20) 



has flown to its IR- fixed point v*. The /3- function asso- 
ciated to the coupling r/ is 



Pv 



d_ 



VR 



VR{~S+{l-d/2)i^— 



InZ + fi—- 



InZA (21) 



Suppose now that rj has a nontrivial fixed point rj* for 
r]ii > 0, i.e. (3ni7]*,b*) = and ^Mm,b*)\^H=r,' > 0- 
(We will show below that at leading order such a fixed 
point exists.) We now express M^l InZ by its value at 
the IR-fixed point bj^ = b* 



2{iy* 






InZ 



(22) 



k J p 



lj-2p-^Qiik+p)(rix,t)-riy,ty) 
,^^kp{C{5xd)+C{5y,t)) 

-subdominant 



We 



(16) 



In order to obtain a long-range term a pole at k+ p — 
is necessary. For d > 2 however, the expression is an- 
alytic. No long-range term is generated. This is easily 
generalized to any order in perturbation theory. 
We now use the fluctuation-dissipation theorem 

e(t - t') (r„ [x, t)rp (y, t')) = (r„ {x, t) (r • V) {y, t')) , 

(17) 

which is derived along the same lines as in [E6[ and which 
we write down in Ito-discretization. (For other discretiza- 
tions, additional terms appear on the r.h.s. which cancel 
the contraction of f and r with the same argument.) It 
is valid for bare and renormalized quantities. Inserting 
the definition of Vafi we conclude that 



Zr, = l 



In the parameterization given above, Z\ is 
Zx = l + 0{7]r) 



(18) 



(19) 



as it has to vanish for 77^^ = 0. These relations are suf- 
ficient to completely solve for the anomalous exponents. 



can then solve for u-^\ h\Z\ 



d^i 



IuZa ^5+{2-d){v* -v) 



(23) 



The exponent z is given by 



z={D + 2v* - ^i— 



InZA)/^^* 



The last two equations can be combined to give 
z = d . 



(24) 



(25) 



This relation is valid as long as 77^ has flown to a non- 
trivial fixed point rj* > 0. We will first study the stability 
of the fixed point rjfj = 0, before analyzing potential fixed 
points for rjB. > 0. Without selfavoidance this is simply 
the line with S — 0, see figure 1. Selfavoidance however 
modifies this phase-separation line. To see this, look at 
Pjj at 7y_R « 0: 



p^^{-S+id-2)iiy*-i^))m + 



4 



(26) 



where we used the fact that fj.-^\nZx = 0{riR). The 



stability condition for the fixed point rjn = is therefore 
6 < {d - 2){i^* - i^) . (27) 



At 1-loop order, the separating line is 
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S={d-2)- 



(28) 



Numerical evaluation yields the thin line separating the 
regions with z = d and z = 2 + D/v* . There is however 
a priori no reason to trust this estimate for membranes, 
i.e. e = A. We know however that in any dimension the 
Flory-estimate t^piory = (2 + D)/{2 + d) is quite a good 
approximation for ly* in the fractal phase, for polymers 
as well as for membranes [ |l2|jl^ ]. Inserting this relation 
we obtain for the separatrix 



d=2(D + l) . 



(29) 



(This is the fat line between the regions with z = d and 
z = 2 + D/v* .) Let us stress that we only use the Flory- 
approximation to estimate v* ^ but not any of the sys- 
tematically wrong assumptions which have to be used to 
derive it. 

Another possibility to get (E9|) is to demand that the 
value of z is continuous on the phase separation line. The 
equivalence of the results obtained by the two methods is 
a consequence of the general structure of the renormal- 
ization group. 

We also can give a rigorous bound for the phase sepa- 
ration line. As i^* < 1, hydrodynamics is always relevant 
for 



d< 



A- D 



(30) 



We still have to check that /3^ has a fixed point for 
ry^ > 0. At 1-loop order Zx is 



Za = 1 



,_.f+0(,J) 



(31) 



The diagram on the r.h.s. is the contraction of the hy- 
drodynamic interaction only. Explicitly this is 



5al3 = 
^ J x<^^^ J k 



Sgf} kakp \ -\xf-D 



= 5a(: 



k^ 

d~ 1 1 



/fc4 



2did-2)6 
The residue is thus positive 



T/^ 



-(5 



5-1 



>0 



(32) 



(33) 



This ensures the stability of the fixed point at least for 
small S. 

In conclusion: We have shown that the dynamical field 
theory (O) for polymerized tethered membranes including 
hydrodynamics is renormalizable and that the dynamical 
scaling exponent z is given to all orders in perturbation 
theory by d. 



It is a pleasure to thank F. David, H. W. Diehl and L. 
Schafer for stimulating discussions. 
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